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Abstract 

The aim of this paper is to understand how to measure the V0 2 and VC0 2 variabilities in indirect calorimetry (IC) 
since we believe they can explain the high variation in the resting energy expenditure (REE) estimation. We 
propose that variabilities should be separately measured from the V0 2 and VC0 2 averages to understand 
technological differences among metabolic monitors when they estimate the REE. To prove this hypothesis the 
mixing chamber (MC) and the breath-by-breath (BbB) techniques measured the V0 2 and VC0 2 averages and their 
variabilities. Variances and power spectrum energies in the 0-0.5 Hertz band were measured to establish technique 
differences in steady and non-steady state. A hybrid calorimeter with both IC techniques studied a population of 
15 volunteers that underwent the clino-orthostatic maneuver in order to produce the two physiological stages. 
The results showed that inter-individual V0 2 and VC0 2 variabilities measured as variances were negligible using the 
MC while variabilities measured as spectral energies using the BbB underwent 71 and 56% (p < 0.05), increase 
respectively. Additionally, the energy analysis showed an unexpected cyclic rhythm at 0.025 Hertz only during the 
orthostatic stage, which is new physiological information, not reported previusly. The V0 2 and VC0 2 inter-individual 
averages increased to 63 and 39% by the MC (p < 0.05) and 32 and 40% using the BbB (p < 0.1), respectively, 
without noticeable statistical differences among techniques. The conclusions are: (a) metabolic monitors should 
simultaneously include the MC and the BbB techniques to correctly interpret the steady or non-steady state 
variabilities effect in the REE estimation, (b) the MC is the appropriate technique to compute averages since it 
behaves as a low-pass filter that minimizes variances, (c) the BbB is the ideal technique to measure the variabilities 
since it can work as a high-pass filter to generate discrete time series able to accomplish spectral analysis, and 
(d) the new physiological information in the V0 2 and VC0 2 variabilities can help to understand why metabolic 
monitors with dissimilar IC techniques give different results in the REE estimation. 

Keywords: V02 and VC02 variabilities; Gas exchange variability; Variability; Open circuit hybrid calorimeter; V02 and 
VC02 power spectrum 



Background 

Indirect calorimetry (IC) has been considered by physicians, 
clinical nutritionists, and researchers as the gold standard 
to estimate the resting energy expenditure (REE) and the 
metabolic substrate utilization in humans (Ferrannini 1988; 
Branson and Joahanningam 2004). The IC technique has 
been implemented through the respiratory gas exchange 
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measurement in order to calculate mainly the averages of 
the oxygen consumption (VO2) and carbon dioxide pro- 
duction (VC0 2 ) using different kinds of metabolic moni- 
tors with dissimilar instrumental techniques that can 
possibly produce different results in the estimation of the 
REE (McClave and Snider 1992). 

The general clinical view assumes that the REE esti- 
mation inconsistency is the result of uncertainties in the 
instruments at the time they compute de V0 2 and VC0 2 
averages in ambulatory or critical care patients (Matarese 
1997). However, this high REE variation is not only con- 
sequence of instrumental precision problem since it has 
different causes other than sensors and electronics noise. 



Springer 



© 201 3 Cadena-Mendez et al.; licensee Springer. This is an open access article distributed under the terms of the Creative 
Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Cadena-Mendez et al. SpringerPlus 2013, 2:688 
http://www.springerplus.eom/content/2/1/688 



Page 2 of 1 2 



We believe that the variation in the REE estimation 
mainly depends on how the VO2 and VCO2 variabilities 
are detected and measured during the standard V0 2 and 
VCO2 average estimations. The variabilities basically 
depend on the patients' physiological state, since each 
one responds differently, during the IC studies, to the 
stimuli produced by their diseases, environmental stress 
conditions (temperature, noise, humidity, etc.) and to 
the metabolic monitoring connector devices such as face 
masks, canopies and mouth connectors (Garby and 
Lammert 1994). 

The problem is that the commercial metabolic moni- 
tors do not separate and measure the VO2 and VCO2 
variabilities from their average computation in order to 
correctly interpret whether these variabilities are noise, 
artifacts, or if they can contain physiological information 
that may help to understand the variation in the REE 
estimation (Karsegard et al. 2010; Sundstrom et al. 2013). 
How to measure the VO2 and VCO2 variabilities in steady 
or non-steady state is an open question since metabolic 
monitors have been designed only to compute the VO2 
and VCO2 averages in IC clinical studies which respond 
to a different type of patient needs such as the pediatric, 
ambulatory, intensive care or those subjects in free move- 
ment conditions (Simonson and DeFronzo 1990). 

The premise is that the metabolic monitors use differ- 
ent instrumental IC techniques which can filter or 
enhance the VO2 and VCO2 variabilities depending on 
the design of the pneumatic hardware. For instance, in 
calorimeters with mass exchange, which can be con- 
sidered as an open-design pneumatic system, the VO2 
and VCO2 measurements of averages are carried out fas- 
ter even in a portable manner but with more sensitivity 
to detect the embedded variabilities (Brown et al. 1984; 
Myers et al. 1990). On the other hand, in metabolic 
monitors with closed-design pneumatic circuit the VO2 
and VCO2 variabilities are filtered because their frequency 
bandwidth response is probably limited by the great size 
of their instrumental set-ups (Holdy 2004). In both cases, 
no study exists that would suggest which type of pneu- 
matic hardware would be more suitable to measure the 
variabilities separately. Thus, metabolic monitor manufac- 
turers are free to promote their own hardware designs, 
which have the sole purpose to obtain the best estimation 
of the V0 2 and VC0 2 averages. 

Currently most of the commercial metabolic monitors 
in use are open circuit systems with the capability to 
perform IC studies in the range of 15-20 minutes. The 
bench performance in their accuracy has been evaluated 
in the range of 2 to 3%, which does not explain the high 
variation in the REE estimation, in critical care and 
ambulatory patients (Reeves et al. 2004 and Cooper 
et al. 2009). The REE variation has been reported as 
inter-subject or intra-subject variability despite following 



the clinical guide for standardizing IC studies. Therefore, 
we assume this guide is only focused on controlling the 
patient's steady state in order to facilitate the prediction 
of the REE in 24 hours without taking into account un- 
steady patients (AARC Clinical Practice Guideline 2004). 
The second assumption is that the IC studies in non- 
steady state, as part of the clinical routine, are accountable 
for the high variability in the REE estimation (Wessel et al. 
1979). Then, current guidance has a limited scope since it 
only describes measurement conditions when patients are 
in steady state such that their physiological reaction to 
thermal and noise environment is monitored and con- 
trolled only through the measurement of the variation 
coefficient (VC) that should be maintained by no more 
than ± 5% in any single 5 minute measurement interval 
(McClave et al. 2003). 

Consequently, any proposal for a new IC clinical guide 
with the capability to define studies in steady and non- 
steady state should consider the following facts: (a) only 
54% of ambulatory patients who are submitted to REE 
estimation comply with the current VC defined thresh- 
old (McClave and Snider 1992; Damask et al. 1983). (b) 
Different IC studies in ambulatory subjects have repor- 
ted intraindividual REE variability of 12.5% and 23% 
during two consecutive days and over a period of 2 years, 
respectively (Mindy et al. 1986). (c) Studies in critically ill 
patients have reported inter-individual REE variability up 
to 64% with different metabolic monitors (Schadewaldt 
et al. 2013). (d) Investigations have shown increments 
beyond 15% in the VC due to what is called physiological 
reactivity which is associated to the effect of the connect- 
ing devices such as face masks, mouth connectors and 
canopies (McAnena et al. 1986; Segal 1987; Isbell et al. 
1991; Forse 1993). (e) Metabolic monitors for patients in a 
hospital or the ambulatory environment having different 
pneumatic techniques and algorithms have not been eval- 
uated to perform IC studies in steady and non-steady state 
conditions, and (f) the most frequently used open circuit 
calorimeters have two different design concepts known as 
the breath by breath (BbB), and the mixing chamber (MC) 
techniques that may produce different REE estimation 
results (Faver et al. 1998; Arch et al. 2006). 

Observations, such as those previously mentioned, 
have triggered our main assumption that the V0 2 and 
VCO2 variability analysis can help to understand the 
variation in the REE estimation. Therefore, we propose 
that variabilities should be separately measured from 
their VO2 and VCO2 averages with the purpose to 
identify the monitors' technological differences in the 
measurement of the variabilities and to search for pos- 
sible physiological information in them during steady 
and non-steady IC studies. To prove this hypothesis, the 
VO2 and VCO2 averages and their variabilities were 
simultaneously measured by using the MC and the BbB 
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techniques. Specifically, variances were used with the 
MC technique and power spectrum functions in the 
0-0.5 Hertz band were constructed and measured in 
the BbB method. The gas exchange in the first 
method was sampled every 20 seconds while the gas 
exchange in the second technique was sampled breath 
by breath in order to generate a stochastic process meas- 
urement as a manner of a discrete time series. Both gas 
exchanges were measured in two consecutive 15 minute 
windows in order to obtain enough data to measure 
trends and to have high frequency resolution when explor- 
ing cyclical rhythms. 

The objective of this work was to develop a pilot study 
where young healthy volunteer population was submitted 
to the physiological clino-orthostatic maneuver (COM) 
with the idea to generate steady/clinostatic and non- 
steady/orthostatic stages as a manner of an instrumental 
bench test (Gonzalez et al. 2013; Cadena et al. 2010). 
A hybrid calorimeter with the MC and the BbB 



techniques was applied to compare inter-individual 
variability changes during the COM stages. The inter- 
individual variabilities' energy measurement was con- 
sidered as appropriate since it can reflect physiological 
phenomena. Then, variance averages {SDv 0 2, SDycoi) 
were used as time energy measurement in the MC tech- 
nique and power spectrum functions (vV02(f), vVC02(f)) 
were used as frequency energy measurement in the BbB 
method. Finally, the inter- individual V0 2 and VC0 2 aver- 
ages were computed with the purpose to compare the 
techniques' performance in the REE estimation. 

Results and discussion 

Variabilities' energy by power spectrum functions 

The Figure 1 shows the BbB inter-subjects VO2 and 
VCO2 variabilities in terms of their averaged power 
spectrum functions vV02(f) and vVC02(f) with units in 
(ml/bth) 2 /Hz in the frequency range of 0-0.5 Hertz. 
Thus, individual spectrum functions were generated 
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Figure 1 The inter-subject variabilities analysis by power spectrums in the 0-0.5 Hertz frequency band. The averaged power spectrum 
functions vV02(f) in graph (a) and the function vVC02(f) in graph (b) show a cyclic rhythm at 0.025 Hz only during the orthostatic stage of 
the COM. 
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when each subject (N = 15) was submitted to the steady/ 
clinostatic and then to the non-steady/orthostatic stages 
of the COM. 

The averaged functions vV02(f) and vVC02(f) showed 
that their main energy was concentrated in the low 
frequency (LF) and medium frequency (MF) bands, 
where the analysis bands were defined as LF = 0-0.04, 
MF = 0.04-0.15 and HF = 0.15-0.50 Hertz. Here, it is 
noteworthy that the band definition for the variabilities 
was analogous to how the heart rate variability is 
analysed with the purpose to facilitate their physiological 
interpretation (Task Force of the European Society of 
Cardiology 1996). The Figure 1 clearly shows that the 
vV02(f) and vVC02(f) have an unexpected cyclic rhythm 
energy with central frequency at 0.025 Hz that is 
generated only during the non-steady/orthostatic stage of 
the COM with an energy increment of 71% and 56%, 
respectively, when this is compared with the steady state 
energy. 

The Figure 2 shows the variabilities' total energy 
(ml/bth) 2 increments by frequency bands. For instance, 
the total energy of vV02(f) was computed as the integral 
below the function so that in orthostatic stage is approxi- 
mately 50% higher in magnitude than the total energy of 
the vVC02(f) for the same stage. A comparative analysis 
for the variability energies as a consequence of the COM 
application is shown in the Table 1. 



56% for the vVC02(f) (p < 0.05) as it is also seen in 
Figure 2. 

Comparative analysis of averages 

The Table 2 shows the VO2 and VCO2 average com- 
parison according to the model of the Figure 3. A 
discrete gas exchange (V02[n] and VC02[n]) model was 
considered in order to compute separately the averages 
for the MC and BbB techniques. Then, the averages for 
the MC technique ( V02[n\ MC and VC02[n] MC ) were 
incremented 63% and 39% (p < 0.05), while the averages 
for the BbB technique (1^02[«] and VC02[«]) were only 
increased 32% and 40% (p < 0.1), when the COM was 
applied. The statistical differences analysis was based on 
the Welch f-test for unequal variances. 

Although all computed averages in Table 2 show 
independent significant statistical differences due to the 
COM, the Figure 4 shows no differences between the IC 
techniques. The post-hoc statistical analysis corroborates 
none differences (p > 0.1) using a Welch f-test for the 
inter-individual V0 2 and VC0 2 averages, after applying 
an ANOVA one-factor test for multiple measurements. 
However, it is worth to observe that the difference 
between the V02[n] MC and the V02[n] averages shows 
a slight tendency (p = 0.39) to reject the null hypothesis 
only during the orthostatic stage. 



Comparative analysis of the variabilities' energy 

The Table 1 shows the comparative analysis of the 
variabilities' energy obtained by the MC and the BbB 
techniques application. The variances averages (SZV02) 2 
and (SDvcoi) 2 are compared against the total varia- 
bilities' averaged energy obtained from the vV02(f) and 
vVC02(f) functions. 

The variance averages increment due to the COM 
were not statistically significant (p > 0.1), while the spec- 
tral total energy increment were 71% for vV02(f) and 



Discussion 

In view of the results, the hypothesis is proved in the 
sense that it is necessary to analyze the VO2 and VCO2 
variabilities in order to understand the causes of varia- 
tions in the REE estimation. The Table 1 shows the 
inter-individual variabilities' energies computed as vari- 
ances which did not show statistical differences during 
the COM test. These unchanged variances lead to the 
interpretation that the MC technique is a better method 
to perform IC studies with patients in steady and non- 
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Figure 2 The inter-subject total energy increments by frequency bands. The increments in the total energy were 71% for the vV02{f) 
56% for the vVC02{f) due to the orthostatic stage. 
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Table 1 Variability energy analyses by technique 



N = 15 subjects 


(SD V02 f MC (ml/bth) 2 


[SD VC02 ) 2 MC (ml/bth) 2 


Total energy BbB (ml/bth) 2 


Total energy BbB (ml/bth) 2 


STEADY/CLINO (SD) 


3.81 (2.6) 


2.96 (2.1) 


35.6 (31.9) 


23.9 (24.5) 


UNSTEADY/ORTHO (SD) 


3.92 (2.5) 


3.37 (4.0) 


60.2 (43.8) 


37.5 (35.2) 


f-test 


p = 0.88 


p = 0.77 


p < 0.05 


p < 0.05 



Variances are for the MC while total energies are for the BbB techniques during the application of the clino-orthostatic maneuver. 
Data in bold means averages and data in parenthesis means standard deviations. 



steady state since the REE is much less affected by the 
VO2 and VCO2 variabilities. The MC technique works 
as a low pass filter such that it suppresses the high 
frequencies of the variabilities generating more adequa- 
cies to measure the averages of the VO2 and VCO2. 
Thus, the analysis of averages in Table 2 is consequent 
with the MC technique performance when an increment 
of 63% for the V02\ 

^\mc ^ s observed while the V02[n] 
only shows an increment of 32%. Additionally, this result 
can be interpreted as the MC technique having the right 
sensitivity to faithfully follow any physiological low fre- 
quency change that affects the V0 2 and VC0 2 averages 
(Bruce 1996). 

These outcomes lead us to understand why the old 
instrument Delta Track II (Datex Finland) has been 
accepted as the reference instrument when new meta- 
bolic monitors are compared against its performance, 
mainly during IC studies in critical care patients. It is 
clear that the canopy in the Delta Track II performs as 
an open circuit MC technique with the capability to re- 
ject high frequency variabilities. Deltra Track's hardware 
can be modelled as a time average filter so that the V0 2 
and VCO2 variabilities from unstable patients produce 
minimum inter-individual variations in the REE estima- 
tion (Severine et al. 2013; Miodownik et al. 2000). 

On the other hand, the power spectrum functions in 
Figure 1 show how IC studies in steady or non-steady 
state can be separated using their energy computation. 
The inter-individual variabilities' energies in Table 2 
show increments of 71% for the vV02{f) and 53% for 
the vVC02(f) when subjects are submitted to the ortho- 
static stage. These results suggest, the BbB technique is 
more suitable for monitoring and controlling the sub- 
jects' physiological condition. This is an improvement 



Table 2 The MC and BbB averages analysis during the 



steady/clinostatic and non-steady/orthostatic stages 



N = 15 subjects 


V02[n] MC 
(ml/min) 


VC02[n] MC 
(ml/min) 


V02[n] 
(ml/min) 


VC02[n] 
(ml/min) 


STEADY/CLINO (SD) 


167 (48) 


1 35 (42) 


152 (37) 


124 (36) 


UNSTEADY/ORTHO (SD) 


273 (74) 


188 (70) 


202 (56.6) 


1 74 (64) 


Welch t-test 


p < 0.05 


p < 0.05 


p<0.1 


p<0.1 



Data in bold means averages and data in parenthesis means standard 
deviations. 



over only using the traditional concept of CV. In 
addition, the cyclical rhythms found at 0.025 Hz in 
Figure 1 can be used as a new figure of merit to describe 
how stable or unstable the patient is during an IC study. 
Although this finding can be used as new physiological 
control information, the interpretation of its origins 
needs more research work. One first approach was 
carried out by dividing the energy of the vV02(f) and 
vVC02(f) in frequency bands similarly to the way that 
heart rate variability (HRV) is processed; after which, 
one second step would be to correlate the energy found 
in LF and MF with the LF energy of the HRV in order to 
discard whether or not the rhythmicity is due to sympa- 
thetic neural control or not (Taylor et al. 2001; Satue 
and Mendez 2012). 

Conclusions 

New and advanced IC metabolic monitors designs 
should consider the following issues: (a) the VO2 and 
VCO2 variabilities should be separately measured 
using the MC and BbB techniques simultaneously in 
order to carry on IC studies in steady or non-steady 
state and to distinguish the origin of the variation of 
the REE estimation, (b) The MC technique is the 
appropriate method to estimate the VO2 and VCO2 
averages, whereas the BbB technique is the most suit- 
able procedure to provide physiological information 
to determine how stable or unstable an IC study is. 
(c) The cyclical rhythm in the BbB technique cannot 
be interpreted as instrumental noise since the O2 and 
CO2 transducers are performing independent mea- 
surements in open circuit calorimeters that allow to 
compute the vV02(f) and vVC02(f) separately, (d) 
The discrete gas exchange modeling of the Figure 3 
contributes to understand why the MC technique is a 
low-pass filter or a moving average measurement sys- 
tem while the BbB technique is a high-pass filter able 
to generate a random process measurement system. 
Therefore, both techniques have different effects in 
the REE estimation, (e) Finally, it is important to no- 
tice that the COM physiological test combined with 
the VO2 and VCO2 variability measurements can be 
considered as an input-output paradigm assessment 
in order to search for the metabolic monitors' 
performance. 
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VC02[n]McA S D vco2 



V02[n] = V02 + variability _ 
VC02[n] = VC02 + variability 



Variabilities 





V02[n]+\/V02[n] 

Discrete-TimW Series 
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Figure 3 The V02 and VC02 variabilities definition. The variabilities are defined in terms of the MC and BbB techniques. The SDV02 and 
SDVC02 are for the mixing chamber (MC) and the vV02[n], vVC02[n] are for the breath by breath (BbB) techniques. The discrete gas exchange is 
defined at the mouth level as V02[n] and VC02[n] such that they include their own variabilities. 



Methodology 

Model for the MC and BbB variabilities and averages 
definition 

The Figure 3 shows a hybrid calorimeter model with the 
MC and the BbB techniques where the V0 2 and VC0 2 
variabilities are defined at the output as SD V o2 and 
SD VC o2 and vV02[n] and vVC02[n], respectively. Like- 
wise, the averages are defined as V02[n\ MC and VC02 
[n] MC for the MC technique and V02[n] and VC02[n] 
for the BbB method. 

This model for the MC and BbB techniques con- 
siders a discrete gas exchange at the input of the 



mouth {V02[n] and VC02[n]) in which the varia- 
bilities are implicitly included before they are separ- 
ately measured. The argument n= 1,2,3... stands for a 
discrete time series that represents the breath by 
breath gas exchange during an IC study. The model 
in Figure 5 explains how the discrete gas exchange is 
formed at the alveolar level. 

Model for the alveolar discrete gas exchange 

The alveolar discrete gas exchange is modelled in 
Figure 5. The assumption is that the continuous gas 
exchange at the alveoli (V02{t) and VC02(t)) is sampled 
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Figure 4 Comparative analysis of averages between techniques. The V02[n] MC and the VC02[n] MC averages correspond to the MC 
technique whereas the V02[n] and VC02[n] averages correspond to the BbB method. The graphs (a) and (b) do not show statistical differences 
between computed averages for both IC techniques. 
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Figure 5 The alveoli gas exchange model. It is assumed that the V02(t) and VC02(t) are sampled by the lung's mechanical ventilation to 
generate the discrete V02[n] and VC02[n] at the mouth level. 



by the lung's mechanical ventilation. Then a discrete gas 
exchange V02[n] and VC02[n] is generated when the 
breath by breath instant flow fit) works as a sampling 
function as in Equations (1) and (2). 

V02{t) x f(t) = V T x {FI0 2 {t)-FE0 2 {t)) (1) 
xf(t) = V02[n] 

VC02(t) x/(f) = V T x FEC0 2 {t) xf(t) 

= VC02[n] (2) 

Where: /(f) is the instant expired flow (L/sec). V T is 
the tidal volume (ml) without BTPS (body, temperature, 
pressure, saturated) to STPD (standard, temperature, 
pressure, dry) volumetric corrections in order to pre- 
serve the simplicity of the model, FI0 2 {t) - FE0 2 (t) is 
the inspired-expired oxygen fraction difference and the 
FEC0 2 (t) is the expired CO2 gas fraction. All gas frac- 
tions are in atmospheric percentages (%). 

The products fit) x V02(t) and fit) x VC02(t) generate 
the continuous sampling for the O2 and CO2 uptake dur- 
ing each expired breath with time duration D lt D 2 , ...D n . 
Hence, individual and different breath-by-breath sample 
volumes are produced as V02[Dn] and VC02[Dn]. These 
volumes are computed as in Equations (3) and (4). The in- 
stantaneous products are done between signals analog to 
digital converter (A/D) at the rate of 10 milliseconds per 
sample in order to avoid numerical integration errors and 
to be according with the sampling Nyquist theorem when 
it is assumed signals with bandwidths below 100 Hz 
(Proakis and Manolakis 1998). 



1 r« 

V02[D n ] = — V T x {FI0 2 {t)-FE0 2 {t)) xf(t) 
Vt Jo 

x g(t-D n )dt 

(3) 

1 f D " 

VC02[D„] = — V T x FEC0 2 (t) x/(t) x g(t-D n )dt 
Vt Jo 

(4) 

The g (t — D n ) are continuous gate functions with the 
same time duration Di, D 2 ,..D n that allow synchro- 
nization to integrate the products between the instant 
flow and the instant gas fractions as it is seen in Figures 6 
and 7. Normalized products fit) x V02{t) and fit) x 
VC02(t) are needed in order to match with the fit) 
peak amplitude the V02(t) and VC02(t) values so 
that Equations (3) and (4) should be divided by 1/V T - 
Figure 6 shows a real example how gas fractions signals 
and the expired instant flow signal are synchronized to 
compute each V02[Dn] and VC02[Dn]. 

BbB discrete time series analysis 

The V02[n] and VC02[n] discrete time series in 
Figure 7 are generated when the computations of each 
V02[Dn] and VC02[Dn] are carried out over the con- 
tinuous signal outputs corresponding to the flow, O2 
and CO2 sensors. The constant time delay of 800 msec 
in Figure 6 is for the synchronization between the in- 
stantaneous flow fit) and the time gas fraction signals 
(FE02(t) and FEC02(t)). This time lag depends on the 
sensors' time response and the delay produced by the 
tubing length which utilizes a flow of 150 ml/min to 
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Figure 6 Gas exchange fractions and instant flow signals. The constant 800 millisecond time delay is considered in order to synchronize the 
computation for each V02[Dn] and VC02[Dn] and the corresponding time series generation. 



sample the sensors. The hybrid calorimeter with the 
open pneumatic circuit is sketched in Figure 9. Each 
V02[Dn] and VC02[Dn] value is placed in a time series 
using the sequence S[n] with mathematical proprieties 
that allow the generation of the BbB discrete time 
series according to Equations (5) and (6). 



V02\n] = V02[D„] x S[n-D„] 



VC02[n] = VC02[D„] x S\n-D n ] 



(5) 
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Figure 7 The V02[n] and VC02[n] a discrete time series examples. Two simultaneous values of V02[Dn] and VC02[Dn] are computed at 
6[n-5] and 6[n-34]. The value of each Dn = 1,2, ..,N is placed at the end of the corresponding expired instant flow. 
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Then, the sequence 8 [n - D n ] is used to geometrically 
place each value of V02[Dn] and VC02[Dn] as a series 
of coefficients at the end of each fit) as it is seen in 
Figure 7. Here, the meaning of Dn is extended as a 
dumb variable {Dn = 1,2,. .n) just to be interpreted as an 
index to generate the BbB discrete time series V02[n] 
and VC02[n]. The Figure 7 shows an example of a 
discrete time series from which the vV02[n] and 
vVC02[n] variabilities are computed. The average values 
(V02[«] and VC02[n]) are calculated from the discrete 
gas exchange as in Equations (7) and (8). 

vV02[n] = V02[n]-V0 2 [n] (7) 

vVC02[n] = VC02[n]-VC0 2 [n] (8) 

Where: 

VC0 2 [n)=^J^ n=1 VC02[n) (10) 

The above averages are computed with approximately 
N = 225 breaths, which are equivalent to a data acquisi- 
tion window of 15 minutes. 

An example of the power spectrum analysis of the 
vV02[n] and vVC02[n] is shown in the Figure 8. A 
linear data interpolation function was used to reformat 
the discrete time series V02[n] and VC02[n]. Then, one 
sample per second was used to resample the reformatted 
discrete time series in order to obtain a frequency 
domain analysis in the range of 0.0 to 0.5 Hz. The 



processing window was selected to capture at least 
15 minutes of data so that a Welch power spectrum 
estimator allowed a maximum resolution of 0.005 Hertz 
using 50% of data overlapping. The frequency band 
analysis was defined in three main regions: low frequen- 
cies (LF = 0-0.04 Hz), medium frequencies (MF = 0.05- 
0.15 Hz) and high frequencies (HF = 0.16-0.5 Hz). These 
band divisions are similar to the heart rate variability 
analysis with special emphasis in the LF and MF bands 
since the HF band is assumed to be related with instant 
flow's frequency (respiratory frequency) activity as it is 
seen in the example of the Figure 8. 

Averages and variabilities in the MC technique 

The measurement of the averages and variabilities using 
the MC technique requires modeling the effect of the 
mixing chamber upon the discrete gas exchange V02[n] 
and VC02[n], having the model in Figure 3 in mind. 
The MC averages should be computed as in Equations 
(11) and (12) using a digital moving average which 
depends on the chamber volume and the number of 
breaths that the chamber storages as a manner of pipe- 
line, prior to obtaining one sample average every 20 sec- 
onds. In our case, the hybrid calorimeter has a chamber 
with a volume of 1.8 Liters so that the number of 
breaths storage in the MC, when the patient's respiratory 
frequency is approximately 15 breaths/min, is approxi- 
mately M = 4 in Equations (11) and (12). And the V02 
[n] MC and VC02[n] MC are computed using the criteria 
of 30 averages to smooth enough the gas exchange. 
Thus, each one of the 30 averages is formed with M 
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Figure 8 Example of a vV02{f) power spectrum in (ml/breath) /Hz. A subject underwent the clino-orthostatic maneuver. Graph (a) shows 
the vV02(fl outlined in black that corresponds to the clino stage. The power spectrum function outlined in gray corresponds to the orthostatic 
stage. Graph (b) shows the power spectra for the continuous instant flow fifj. 
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breaths to obtain one average sample during a total 
of 15 minutes per each IC study. 



V02\n] 



V02\n] 



1 

30 



E30 f 1 -c-^M „ . -] 

<^ > V0 2 \n-m\\ 

E30 f 1 
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920 sec 
(11) 



1 ^-^30 

30 



(12) 



The measurement of the variabilities in the MC tech- 
nique was computed as in Equations (13) and (14). Even, 
these equations allow the calculation of the VCs according 
to the clinical practice guidelines as it is shown in (15). 



SD V 



1 / ■ , . ( 1 v-^M 



V02\n-m\ 



(13) 
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Hybrid indirect calorimeter hardware 

A specific open-circuit hybrid indirect calorimeter 
(MGM-3) was designed and manufactured for the pur- 
pose of this work which was based in the design of 
Westenskow et al. (1984). The MC and the BbB 



techniques were fused in the MGM-3 as it can be 
seen in Figure 9. The patient's half mask works either 
by passing the expired gas through the 1.8 L mixing 
chamber to implement the MC technique or by directly 
connecting the expired gas to a hot-wire flowmeter 
(TSI Inc, USA) to implement the BbB technique. 

The MGM-3 calibration and quality control unit was a 
microprocessor based design and was calibrated every 
5 minutes using a reference gas cylinder with a certified 
mixture of 21% 0 2 , 10% C0 2 , complemented with N 2 . 
Additionally, two more gas certified mixture cylinders 
(15% 0 2 , 4% C0 2 and 18% 0 2 , 3% C0 2 , Praxair) were 
used to adjust the transducer offsets and gains for the 
case of doing IC studies in ambulatory patients. The MC 
technique was implemented by displaying values of Ve 
(expired volume minute in L • min" 1 ), RF (respiratory 
frequency in breaths • min" 1 ), V0 2 (ml • min' 1 ), VC0 2 
(ml • min" 1 ), Vt (tidal volume in ml • breath" 1 ) and RQ 
(respiratory quotient V0 2 /VC0 2 ) every 20 seconds. 
These readings were automatically corrected and dis- 
played at STPD conditions after measuring volumes and 
fractions at BTPS conditions (2400 meters above the sea 
level at Mexico City, 590 ± 3 mmHg, and expired gases' 
temperature). The V0 2 was computed using the Haldane 
correction. 

Experimental design and data processing 

A population of 15 young normal volunteer subjects 
without a history of any chronic disease was studied. 
The ages ranged from 18 to 30 years with a body mass 
index (BMI) average of 24.2 ± 3.8 Kg • m" 2 . All subjects 
gave signed informed consent to be studied in the 
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block diagram of the hybrid calorimeter MGM-3 is shown. The MC section is outlined in dotted lines and the blue blocks point 
pneumatic open circuit section. The PC-104 is a dedicated computer to obtain real time data from the flow meter, the 0 2 and 
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morning after 8 hours of fasting. Subjects were asked 
to perform the active clino-orthostatic maneuver 
(COM). First, a 5 minute period of relaxing was used 
before he/she lied down on a couch and was then 
submitted to the COM while connected to the 
MGM-3 calorimeter. Two 30 minutes periods were 
used to implement the measurement protocol: 15 minutes 
for the MC technique and 15 min for the BbB tech- 
nique in each COM position. All of the measure- 
ments were made in the same room maintaining 
constant temperature and data collection by the same 
expert team in all cases. 

Comparative statistical paired data analysis was applied 
intra-groups. The MC averages V02[n] MC and VC02[n] M c 
were compared against the BbB averages V02[n] and 
VC02[n] . Similarly, variances averages (SDyoa) 2 an d 
(SDycaz) 2 were compared against total spectral energy 
averages. The statistics analysis was parametric since the 
variables were considered to be Gaussians, once they were 
tested by the Kolmogorov-Smirnov test. Then, two-tailed 
paired t-tests (Welch version) were used as appropriate for 
unequal variances. In all cases, the null hypothesis was 
rejected when p < 0.1 since this experiment was considered 
to be a pilot study. 

Abbreviations 

IC: Indirect calorimetry; COM: Clino-orthostic maneuver; REE: Resting energy 
expenditure (Kcal/day); V0 2 : Oxygen consumption (ml/min); VC0 2 : Dioxide 
production (ml/min); MC: Mixing chamber technique; BbB: Breath-by-breath 
technique; VC: Variation coefficient (SD/average); SD: Standard deviation; 
SDvoi- V0 2 variance (energy) (ml/breath) 2 ; SD^coi- VC0 2 variance (energy) 
(ml/breath) 2 ; vV02(f): V0 2 variability average power spectrum function 
(ml/breath) 2 /Hz; vVC02(f): VC0 2 variability average power spectrum function 
(ml/breath) 2 /Hz; Total energy: = \\Z.vV02(f)df (ml/breath) 2 ; V02[n] MC : 
V0 2 average by the MC technique (ml/min); VC02[n] MC : VC0 2 average by 
the MC technique (ml/min); l/02[n]: V0 2 average by the BbB technique 
(ml/min); VC02[n]: VC0 2 average by the BbB technique (ml/min); V02[n] 
and VC02[n\. Discrete gas exchange at the mouth; V02(t) and VC02 
(t): Continuous gas exchange at the alveoli; vV02[n]: V0 2 time discrete series 
variability (ml/min); vVC02[n]: VC0 2 time discrete series variability (ml/min). 
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